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abstract 



r\ \ Baryons of lower spin in partially quenched large- A^c QCD are studied with particular emphasis 
cd ' to interpolation between standard unquenched and fully quenched limits. In large Nc limit of 
partially quenched QCD, we calculate Am, the chiral one-loop correction to baryon masses. 
We find that leading order contribution to Am is independent of number of ghost quarks 
introduced. For finite N^ Am does satisfy the Bernard-Golterman's third theorem and has no 
infared quenching singularity except for fully quenched limit. At large Nc limit, however, we 
show that the third theorem is bypassed and non-trivial quenched chiral corrections do arise. In 
unquenched limit, we also show that standard chiral perturbation theory results are reproduced 
in which rj' loop contributions are explicitly taken into account. 



1 Work supported in part by NSF Grant, NSF-KOSEF Bilateral Grant, KOSEF SRC-Program, Ministry of 
Education Grant BSRI 97-2418, SNU Research Fund, and the Korea Foundation for Advanced Studies Faculty 
Fellowship. 



1 Introduction 

With impressive advance in numerical simulation of QCD on a lattice, the goal of first-principle 
calculation for hadron spectrum as well as other physical quantities is achievable in the near 
future. At present, however, nearly all QCD simulations are based on various degrees of approx- 
imation. Therefore, understanding error and deviation incurred by such approximations to the 
calculated physical quantity is of prime theoretical importance. Among the approximations, 
adopted most frequently is the quenched approximation [|l], H. By quenched approximation, 
which has been introduced in order to accelerate generation of the gauge field configurations, 
one drops out quark determinants and keeps only the valence quarks. While the approximation 
is well suited when quark masses are taken sufficiently heavy, extrapolation of quark masses to 
chiral limit might lead to potentially significant effects. In fact, Sharpe P, Q and Bernard and 
Golterman |Q |^ have pointed out that the quenched approximation leads to sizable errors in 
the chiral limit. If this is the case for a generic physical quantity, then one needs to understand 
better the errors incurred by the quenched approximation before a prediction or a conclusion 
is drawn. 

To investigate systematically an error introduced by the quenched approximation in the 
chiral limit, Sharpe p, ^ and Bernard and Golterman [^ have developed chiral perturbation 
theory for quenched QCD. This so-called quenched chiral perturbation theory (QxPT) utilizes 
the trick by Morel of introducing ghost quarks, which have the same masses as the original 
quarks but carries opposite statistics. Thus, the chiral fiavor symmetry is extended to graded, 
super-fiavor symmetry and it is this symmetry which organizes the QxPT. Moreover, diagrams 
with quark loops are cancelled by the same diagrams but with ghost quark loops. Since these 
pioneering works, using the QxPT, chiral behavior of quenched QCD has been studied exten- 
sively for pseudo-scalar mesons 0, ||, |^, baryons [|], vector and tensor mesons [|, |10|, heavy 
mesons [[iT], |l^ , weak matrix elements Bk 0, |13[ , baryon axial charge |]14| , heavy baryons |l5l 



and pion scattering length [0. Numerous lattice data have been compared to results from 
QxPT: see, reviews by Sharpe ^^, Gottlieb |]18[ and Okawa [^ and references therein. 



Through these extensive study, it has been concluded that main effect of quenched approx- 
imation is due to peculiarity of super-r^' meson (direct counterpart of "anomalous" rj'). This 
super-?]' meson has the same mass as the other Goldstone bosons (same single-pole position 
of propagators) and has to be retained in the QxPT, unlike the standard xPT. Peculiar to 
super-?]' is the extra double-pole term present in the propagator, which does not permit a par- 
ticle interpretation at all. Most notably, it has been found that the double-pole term in the 
super-?]' propagator leads to new infrared singular, non-analytic chiral corrections that have 



no counterpart in the full QCD. In order to understand better such pathological behavior, it 
should be desirable to be able to interpolate, if possible at all, QCD between unquenched and 
fully quenched limits and examine changes of physics in the chiral limit. Indeed, with such 
motivation, Bernard and Golterman have developed partially quenched QCD, in which only 
a subset of quarks are paired with ghost quarks [|. By formulating partially quenched chiral 
perturbation theory (PQ^PT), Bernard and Golterman have shown that change of chiral be- 
havior for Goldstone mesons can be indeed understood from interpolation between standard 
QCD and fully quenched QCD. 

In this paper, we study chiral dynamics of large Nc baryons in partially quenched QCD. 
Previously, in the fully quenched QCD, baryons have been studied by Labrenz and Sharpe p|. 
They have studied quenching effect to the behavior of baryon mass spectrum in the chiral limit 
and have observed, most notably, that pattern of the correction is substantially different from 
that in the unquenched case. The leading order chiral correction in QxPT scales like m^^, which 
is more singular than the standard m^/^ xPT correction. Prompted by such interesting results, 
in this paper, by formulating partially quenched chiral perturbation theory for baryons, we 
examine change in chiral behavior of baryon mass spectrum as the QCD interpolates between 
standard and fully quenched limits. In our investigation, we also invoke large Nc limit and 
associated planar symmetry. Combined application of chiral perturbation theory and large Nc 
expansion are expected to constrain the low-energy interactions of baryons with the Goldstone 
bosons more effectively than either method alone. We have found that the large Nc limit 
affects the chiral dynamics in two interesting ways. First, the positions (on the complex plane) 
of the poles of the rj' propagator depend explicitly on the value of Nc- In PQxPT, the rj' meson 
propagator has two simple poles: pion pole at p"^ = m? and a shifted pole at p^ = m? + Mq. The 
"pole shift" Mq originates from the hairpin diagram, which is 1/Nc suppressed. As Nc — * cxd, 
the shifted pole will merge to the pion pole. As we will show below, there arises an interesting 
cancellation between contributions from each of these two poles in the rj' propagator. Second, 
the size of Nc also controls the form of axial current couplings in the baryon sector. The rj' 
coupling to a single quark line inside a baryon is related to its pion counterpart, as vr and rj' 
are related by the "planar symmetry", which is exact in the large Nc limit. This symmetry, 
however, is not manifest at the hadron level. While pion coupling to different quark lines 
interfere constructively, the rj' coupling has huge cancellation. As a result, chiral loops involving 
TT and rj' appear at different order in large Nc expansion. 

This paper is organized as follows. In section 2, after recapitulation of baryons in the large- 



^As Bernard and Golterman point out, partially quenched QCD arises also quite naturally for description of 
lattice QCD with staggered fermions or of lattice QCD with Wilson valence quarks and staggered sea-quarks. 



Nc limit we formulate partially quenched cliiral perturbation theory for large Nf. baryons. We 
pay particular attention to the mass corrections of nucleons and A in the large N^ and chiral 
limits. In section 3, we calculate infrared singular, non- analytic chiral corrections to the baryon 
masses. Moreover, we provide anatomy of chiral one-loop corrections and identify contributions 
of rj' at leading order in 1/Nc expansion. In section 4, we investigate the effect of next leading 
order corrections in 1/Nc expansions, followed by some discussions. 



2 Baryons in Partially Quenched Large- A^c QCD 

2.1 Baryon Dynamics in Large Nc QCD 

Large Nc meson dynamics has been first studied in Ref. ||20[. Baryon dynamics in large- A'"c QCD 



has been studied originally in Ref. [^ and, more recently, in Refs. p^, ^ 0, |^, ^. We 
recapitulate essential aspects of their results that will become relevant for baryons in partially 
quenched QCD (PQ-QCD). For simplicity, we will mainly focus on the case with n = 2 [u 
and d). Since the Nc quarks inside a baryon should form a color singlet, which is completely 
antisymmetric, the spin-flavor part of the wave function must be completely symmetric. As 
a result, the lowest lying baryons have I = J = ^, §,..., which are usually identified as the 
observed states N, A, . . . ^ 

The interaction of baryons with Goldstone bosons are simplified in the large Nc limit as 
the dynamical symmetries are enlarged both in the meson and the baryon sector. In the 
Goldstone boson sector, there arises the planar symmetry (also called "nonet symmetry" in 
the literature), decreeing that the rj' meson should be combined with the other Goldstone 
bosons into a n x n representation of U{n) flavor symmetry group, where n is the number of 



light flavors pO|, ^. More explicitly, for n = 2, Goldstone boson fields are represented by the 

2x2 matrix: 

1 / jLjf. Jil 7r+ \ 

= vr'^T'^ + V-= = 0"T"=h^ _v^ ^_^], (1) 

where T" (a = 1, 2, 3) are the SU(2) generators, 1 is the identity, and T" {a = 0, 1, 2, 3) 
are the U(2) generators, i.e., {T"} = {T"^} U {1}. As we will see below, the Nc ^ oo planar 
symmetry implies that the same set of coupling constants will control both the vr and the rj' 
interactions. 



^See Ref. [M for detailed discussion on validity of such an identification. 



In the large A''^ limit, the Goldstone boson kinetic term is 

C^=^-TT{d^^ad,a) + ..., (2) 

o 

where 

a = cxp{2z(j)/f), (3) 

and ellipses denote higher derivative interactions that are irrelevant for low-energy dynamics. 
The degeneracy between vr and r]' persist even if one endows a small mass to the light quarks 
and breaks the chiral symmetry. The quark mass leads to a non-zero mass ml. = rni, = m?. 
The U(1)a anomaly, which breaks the planar symmetry and renders a heavy mass to r^', shows 
up as an 1/Nc suppressed correction: 

£i = i(m2(Tr0)2 + Ao(9^Tr0)2) =1(^2^^ + A,d,r,'d^i), (4) 

where n is the number of light flavors, and Aq and rn^ are of order 0{N^^). In the standard, 
unquenched xPT, the effect of Ci can be resummed: 



(5) 



Y.— 2\ [nAop'^ + nmQj — 

1 1 

1 — uAqp'^ — {rn^ + n wig/(l — n Aq)) 



One finds that new position of single pole at uii, = m^ + nm'^/{l — uAq). Thus, in the chiral 
limit, m^i ~ 0{N~^), as expected. 

On the other hand, the baryon sector exhibits an SU(2n) spin-flavor symmetry . For 
simplicity, we will restrict ourselves with 2 light flavors for foregoing discussions, in which case 
the spin-flavor SU(4) symmetry will be generated by: 

(6) 

Here, the spin and isospin operators on individual quark lines are denoted by S* {i = x, y and 
z are three spacelike directions perpendicular to the baryon velocity) and T" respectively, and 
1 is the identity operator in spin and isospin spaces. Note that there are 3 J's, 3 I's and 9 



k=i 


r 


fc=i 


Q^a. 


k=i 



Hn the literature, m§ and ^o are often written as n^/3 and a/3 so the for n = 3, mL = m^ + ^^/(l — a). 
In this paper, however, we naturally opted for the notation without the factors of 1/3 as we are mainly working 
with just two light flavors. 

^Our description of the spin- flavor symmetry for large Nc baryons follow closely Ref. [^3[ |2^, ^ . Part of the 
reason is that the formalism in Ref. ^^, which has been widely employed in the literature, assumes unitarity, 
hence, cannot be directly applicable to non-unitary, (partially) quenched QCD we study presently. However, 
the physical predictions of any of these formalisms should be identical. 



G's, making up the correct number of generators Eq. (15) for the spin-flavor SU(4) algebra. 
In what follows, we will adopt a useful notation 

G^" = J2ql{S^0T'^)qk, (7) 

k=l 

which is equal to J' and G*° for zero and non-zero a respectively. 

The axial current couplings of large Nc baryons are simplified by the fact that multi-quark 
operators are suppressed by powers of l/N^. In the leading order of the large Nc expansion, 
the axial currents couple through single-quark operators. Hence, interaction of baryons to 
Goldstone bosons can be expressed as [^ ^: 



/ 



%-. I 5V" Y. ll (S* ® T'^) Ik + dW E 4 (S' ® 1) 9fc + ■ ■ • , (8) 



fc=l fc=l / 

where f = fn = frj' denotes Goldstone boson decay constant, and ellipses denote higher-order 
interactions involving more than one Goldstone bosons. Note that the vr and r]' couplings to 
baryons are taken by the same coupling constant g, as decreed by planar symmetry in the 
large- A'"c limit. 

One can read off from Eq. (^ the meson-baryon-baryon couplings by summing up contri- 
butions from each individual quark lines. The t]' meson couples to individual quark through 
S* ® 1. Hence, rj' meson coupling to a baryon obtained by summing over all constituent quarks 
is given by spin J* of the baryon. That is, 

Nc 

c,,BB = ^f (av) E ^I(s^ ® 1) 9^^ = 4(^V)J^ • (9) 

J k=l J 

Since we are interested only in states with I = J ~ (9(A^°), we will find that the rj'BB coupling 
is of order 0{N~^^'^) (as / ~ (9(A^y^)). Note that C^'be contains sum over Nc quark terms, 
each of them being of order 1// = O^N^^^"^). Hence, one would naively suspect that the sum 
to be of order {Nc/ f) ~ 0{Nc^'^). This reasoning is fallacious, however, as for lower spin states, 
there exists huge cancellations among the individual quark spins. 

One the other hand, pion coupling is proportional to S* (g) T'', and acts coherently over the 



individual quarks |23, 24, 25, EE U%. That this is so can be seen clearly using the operator 



identity |2|, |7 



E E?I(s^®i)?0 +E E^Ki^n^O +4E E^Us^^n^. 

i \k=l / a \k=l / i,a \k=l / 



= J2 + I^ + 4G2 = (|iV2 + 3iV,) 1, (10) 

which is nothing but the Casimir identity for the spin-flavor SU(4). Since both J^ and I^ are 
of order 0{N^), the operator G'"" = J2k=i ll (S' ® T") Qk must be of order 0{N^) so that the 
sum is of order 0{N^). As a result, rj' coupling is one order in l/N^ smaller than the pion 
counterpart, and we can see immediately that r]' loops will be suppressed by l/N^ with respect 
to pion loops. Since all infared quenched singularities arise from r]' loops, we expect these 
infared quenched singularities are subleading in orders of 1/A^c- Below, we will show explicitly 
that this is indeed the case. 



2.2 PQxPT for Large-iVc Baryons 

Now we are ready to (partially) quench the theory. Following the example of Ref. |0, one 
includes in the theory, in addition to the n quarks, k ghosts (spin-| objects with bosonic 
statistics), where k < n. The nxn Goldstone boson matrix (p gets enlarged to an (n+k) x (n+k) 
supermatrix $, 

*-(' i)^ (") 

where x is the nx k matrix of quark-antighost fermionic "mesons", its charge conjugate x^ the 
k X n matrix of ghost- antiquark "mesons" , and the k x k matrix of ghost- ant ighost bosonic 
mesons. The field $ transforms a.s a. {n + k) x [n + k) representation of the enlarged symmetry 
algebra U(n + k) (or more exactly, the l]{n\k) graded algebra). We will denote the generators 
of this \]{n + k) by T^ and 

$ = $^T^. (12) 

Note that {T"} C {T"} C {T^}. 

The lagrangians Cq and £i are generalized to 

Co = ^Strid'^Ed.E) + ..., (13) 

o 

A = ^(m2(Str$)2 + Ao(5^Str<l>)2) (14) 

where 

S = exp(2^<l>//), (15) 

and "Str" denotes the supertrace, the sum of the first n entries on the diagonal subtracted by 
the sum of the last k entries. The propagator for flavored states in $ has a simple pole at 

6 



p^ = m? , while for the flavor- neutral states (states on the diagonal of $) the propagator is 



Gij 



Sue 



ij^i 



p^ — m? An {1 + /\n Aq) p^ — {rn^ + An iTio) p'^ — m?' 



(16) 



The grading index ej is such that e(gj) = +l,e(gj) = —1. Note that the second term (which 
proportional to 1/ An) mixes quark-antiquark meson with ghost-antighost mesons. Both the 
first and the second terms have a simple pole at p^ = m?, but the second term has an additional 
pole, which we will call the "shifted pole". Also note that, as ra — > fc, the shifted pole moves 
back to the pion pole. The propagator will then have a double pole at m?, which is a well-known 
result in QxPT §. 

The partial quenching of the baryon chiral lagrangian is straightforward. The more general 
PQxPT lagrangian in the large N^. limit is 

C = i^(g (a*$^) G'^ + h (S'Str*) 3\ ) . (17) 

The first term is a simple generalization of the unquenched lagrangian, with 

G^^ = E4(S*®T^)?fc' (18) 

k=l 

and g an undetermined coupling constant. The second term describes the "hairpin coupling" 
which couples flavor- neutral states in $ to the baryon operator J^, defined by 

J: = E4(S^®1^)?/^' (19) 

fc=i 

where 1^ is the identity of U(n| A;), with n I's and k — I's on the diagonal. In this paper, however, 
we will set the hairpin coupling constant h to zero to simplify the physics. A discussion of the 
possible effects of a non-zero h will be given wherever appropriate in later sections. 



3 Non-analytic Chiral Correction to Baryon Masses 

3.1 Chiral One-loop Correction 

We are now ready to calculate the chiral one-loop correction to the baryon masses in PQxPT. 
While our theory has n quarks and k ghosts, all assumed to be degenerate, we will focus on 
baryons with just uoi d quarks, though the non (m, d) quarks (which will be collectively referred 



Fig. 1 



Figure 1: Chiral one- loop diagrams to baryon two-point function. Horizontal line denotes a 
baryon and upper internal line denotes Goldstone meson multiplets. 



as strange quarks) may appear in loops. All the diagrams involved have the form as shown in 
Fig. 1 and yield the same one-loop Feynman integral: 

^i('^') = -T^^' (20) 

as those in the standard xPT. We will restrict, in this paper, to the choice of hairpin coupling 
h = 0, the ghost-antighost meson does not contribute, and the only relevant coupling is the 
Lagrangian Eq. (^. 

Depending on whether the Goldstone boson couples to the same quark line at the two 
vertices or not, there are two different types of contribution to the mass correction Am. 

(i) 1 quark line chiral correction: this class of corrections arises when in which the 
Goldstone meson couples to the same quark line at both vertices. The quark level diagrams are 
Fig. 2a and 2b. Notice that Fig. 2b coupling is possible only when the intermediate Goldstone 
meson state is rj'. The contribution of this class is given by 



Am = g' Yl [ Mm')-{ql{S'®T^)qk 



k=l 



2 



+ ^•(Xi(M2)-Xi(m2))-(gi(S^®l)g,)']. (21) 

The first line comes from the chiral loops with flavored Goldstone mesons, and the flavor-neutral 
mesons propagating with the first term in propagator Eq. ([T^). Note that the sum over \]{n\k) 
flavor generators T^ originates from the internal loop in Fig. 2a, which can be either a quark 
or a ghost. The second term line originates from the second term in Eq. (|T6|) , with 

, 9x 1 fni? + An-m'^\ , , 

Xi M' = ■ Xi ——- -^ . (22 

8 



Fig. 2a 



Fig. 2b 





Fig. 2c 



Fig. 2d 



Figure 2: Quark line diagarm for chiral one-loop correction to baryon two-point function, (a) 
Goldstone boson coupling to the same quark line , (b) hairpin propagator to the same quark 
line, (c) Goldstone boson coupling to the different quark lines, (d) hairpin propagator to the 
different quark lines. 



(ii) 2 quark line chiral correction: this is the class of corrections in which the Gold- 
stone meson couples to different quark lines at each vertices. The quark level diagrams are 
depicted in Fig. 2c and 2d. The contribution is given by 

Am = / E'[ Um') ■ {qi (S^ ® T") q,) ■ (g] (S^ ® T") g,) 
k,e=i 

+ (^)-(Ji(M2)-Xi(m2))-(gi,(S^®l)gfc)-(gl(S*®l)g£)]. (23) 

The prime on the summation sign denotes the condition k ^ i. Interpretation of different terms 
is identical with that of (i). Note that in the first line, even though the sum should be in 
principle over all possible Goldstone bosons in $, there is actually no internal quark/ghost loops 
in Fig. 2c, hence, the ghosts do not contribute. As a result, we can just sum over intermediate 
states in (p instead, leading to a sum over T" instead of T"^. 

Now it is time to perform the sums over the quark lines. First, we have terms involving the 
SU(2) generators T*^, i.e., the pion loop contributions. 

Am^_ioop = /JiK). 5] (gt(S^®T")gfc)(gl(S^®T'^)g,) 



kl=l 



where we have used the identity p5 



Nc 



Nc 



Y: {qi (S^ ® T") q,) [qj (S^ ® T») qi) = (Yl i (S* ® T») q, 

k,l=l 



(24) 



|iV2 + fiV, - ij2. (25) 



\k=l 



Second, we have terms involving the flavor identity 1, i.e., the t]' loop contributions. 



Am„' 



»?' — loop 



/\n. ^ 



9 



An ^ 



Nc 
kl=l 



Y (gI(S^®l)g.)(gl(S^®l)g.; 
J', (26) 



where we have used another identity, again from p5 

N, 



Y: (gl (S^ ® 1) gfc) (gj (S^ ® 1) g,) = E ?I(S^ ® 1) ^0 = ^' 



k,l=l 



\k=l 



(27) 



Finally, there are terms involving T"^, those T"^ which do not belongs to the set of T". 
Physically these come from Fig. 2a, with a ghost or a strange quark running around the internal 
loop. There are n — 2 strange quarks and k ghosts, so the contribution is 



ZA771/^^-^_1qqP ylATl ^j^'^X— loop with 1 strange quark- 

10 



(28) 



The combinatorics factor of the latter has been calculated, in Ref. [22, 25] to be |A^c- Hence, 



Amx/x-ioop = g'Ti{m'){/\n - 2){N,. (29) 

Combining all the contributions, we have chiral correction to baryons in partially quenched 
large N^ QCD: 

Am = + /[ (|iV2 + |AniV, + ij2).Ji(m2) + 2-^-j2.(ji(Af2)-Ji(m2)) ]. (30) 

Note that the final result depends only on An, which counts the difference in the number of 
physical quarks n and the number of ghost quarks /c, but not on n and k seperately. Physically 
this reflects the fact that physical quantities should not be changed upon introduction of an 
extra set of degenerate quark and ghost quark pair, as the contribution should cancel out 
completely. To sum up, Eq. (^) is the mass correction to large N^. baryons in PQxPT. Note 
that this single expression works for all states with spin J as long as J <C Nc- 

Before we move on, we note that our result can be expressed in terms of just Ti, the same 
functional form as standard xPT results. The "new chiral singularities" or "quenched infared 
divergences" in QxPT do not appear except the limit An — > 0. This is in agreement with the 
Bernard-Golterman's third theorem , which states that quenched infared divergences appear 
if and only if one or more of the valence quarks are fully quenched. Since the theory we are 
considering is only partially quenched, we do not see the new chiral singularities. 

3.2 The Large-A^c Decomposition of Chiral Corrections 

So far, we have analyzed the chiral one-loop correction to baryon mass spectrum without explicit 
reference to large- iVc counting. In this section, we will decompose the chiral one-loop correction 
explicitly as a function of three distinct parameters: chiral symmetry breaking parameter m^, 
planar symmetry parameter l/N^ and quenching parameter An. 

The chiral mass correction Eq. (^) contains three terms with different N^ and An depen- 
dences: 

Am = Am(+i,o) + Am(o,+i) + Am(_i,o) + Am(„i „i), (31) 

with 



Am(+i,o) =l-9'-N^- Mm') ~ O [N^ An') , (32) 

Am(o,+i) =l-g'-An-N,-I,{m^) r^O(N^,An'), (33) 
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Am(_i,o) 



Am 



(-1-1) 



= \-g^i^-X,{m') 




~C(Ar-\An°) 


(34) 


g^3^-£.{UM^)- 


'Um')) 


r^ O (^N~\ An-^) . 


(35) 



where we have pulled out implicit l/Nc dependence through l/p in our definition of Xi. We 
immediately observe that the terms that depend on J^, which contribute to A-N mass splitting, 
appear only at order 0{N~^). This is in accordance with the large Nc counting given in 
Ref. I2g,||,|2|,g5|,g6|,0. 



We can get a better understanding by identifying the diagrams behind these four terms. The 
first term Ar7i(_|_io), is the dominant contribution in the large Nc limit. It comes from Fig. 2c, 
where a pion attaches to different quark lines at the two vertices. There is a combinatoric 
factor Nc at each vertex. As a result the sum is proportional to N"^/ p ~ A^^- Since there is 
no quark loop in these diagrams, they are completely unaffected by the introduction of ghosts. 
Hence we come to the conclusion that, at leading order of 1/Nc, chiral one-loop corrections 
to baryon masses are independent of the degree of quenching. Fig. 2c contributes at both the 
fully quenched and the unquenched limits, in sharp constrast to the PQxPT for mesons, where 
the diagrams appear in the fully quenched limit always vanish in the unquenched limit, and 
vice versa. In passing, we note that the higher order correction Ar7i(_i^o) also comes from the 
Fig. 2c. 

The term Am(o,+i) comes from Fig. 2a, where the pion attaches to the same quark line at 
both vertices. As a result, the combinatoric factor is Nc, and the whole term is of the order 
Nc/ P ~ N^. Fig. 2a has an internal loop, which may be a quark or a ghost. As a result, it is 
sensitive to quenching, hence, proportional to An. Lastly, Ar7i(_i _i) is the rj' loop contribution. 
It is l/A'c suppressed with respect to the leading pion loop correction as discussed before, and 
the 1/An factor comes from the hairpin term in the propagator. Actually, Ar7i(_i _i) contains 
imphcit suppression factors as the hairpin diagrams are OZI, hence, l/Nc suppressed. We will 
study this in detail in section 3.5. 

3.3 Quenching-Senstivity of Leading Order Chiral Corrections 

In the previous section, we have shown that at leading order of l/N^ chiral one-loop corrections 
to baryon masses are independent of An. A casual reader may suggest that this result is a 
rather trivial observation. Isn't it true that, since internal quark loops are suppressed by l/Nc 
pO| , the effect of quenching must vanish in the large Nc limit uniformly for any degree of 



quenching, hence, our result trivially follows? 
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The above line of reasoning actually turns out fallacious. While it is true that quark loops 
are suppressed by l/A^c? it is definitely not true that the chiral corrections are independent of 
the degree of quenching at leading order correction of l/N^. expansion. There are numerous 
counter-examples that exhibits large Nc corrections that depend sensitively on the degree of 
quenching. As mentioned above, the mass of the p meson receives a chiral correction which is 
proportional to m? in unquenched chiral perturbation theory, but this piece is absent in the 
quenched theory H], and is in fact proportional to An in the partially quenched theory |TD 



Similarly, the chiral correction to heavy meson mass is non-zero in the unquenched theory but 
vanish in the fully quenched limit [[ill , |12| . 



To see flaws in the above argument, note that chiral corrections always come from one 
Goldstone boson loop level. For the case of the p meson and heavy mesons, all Goldstone 
boson loop diagrams always contain a quark loop. Hence the quark loop enter at leading 
order of the l/Nf. expansion, and hence the chiral corrections are sensitive to the degree of 
quenching. On the other hand, there are Goldstone boson loop diagrams contributing the 
baryon mass correction which does not contain any internal quark loop (Fig. 2c), which is also 
the leading order contribution in large Nf.. As a result, the leading order result Am(+i q) is 
An-independent. 

From the above discussion, it should be clear that the question whether the leading order 
(in 1/Nc) chiral correction of a given matter field is quenching-sensitive is highly non-trivial. It 
simply is not true that for the chiral corrections for mesons are always quenching-sensitive while 
those for baryons are quenching-independent. In fact one can show that the chiral corrections 
of flavor non-singlet mesons and baryons with two or more light quarks are always quenching- 
sensitive [p8 |. 



Lastly, we will discuss the implication of the quenching-independence of Am(+i,o)- We have 
shown that the effects of quenching are 1/A'"c suppressed. While we cannot prove that these 
1/A^c suppressed effects are actually negligible (they may come with huge coefficients), it is a 
well-known rule of thumb in hadron phenomenology that 1/A^c corrections are usually small in 
comparison with the leading order result if the latter is non-vanishing. In this case, this rule of 
thumb provides circumstantial evidence that the chiral corrections to baryon masses are small, 
and gives us more confidence in the utility of quenched QCD as an approximation of the real 
world of QCD. 
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3.4 Fully Quenched QCD Limit: An = 

As recalled earlier, the partially quenched QCD can interpolate and bridge between the two 
extreme limits, viz. fully quenched and unquenched QCD theories by dialing the analytic 
parameter An between zero and n. For example, one would expect that by setting An = the 
previously known QxPT results should be recovered. In particular, strongly infared singular 
non-analytic corrections should reappear. In this section, we will study the fully quenched limit 
carefully to understand how these quenched chiral singularities arise. 

First of all, we note that, in the fully quenched limit. An -^ 0, the contributions Am(+i q) 
and Am(_i 0) remain intact and A?7i(o,+i) vanishes identically. On the other hand, the con- 
tribution Ani(„i _i) apparently diverges in this limit. Note, however, that An7.(„i _i) can be 
re-expressed as 

2 2 2// 1 \ T /"^^ + An-mQ^ ^ 



As An -^ 0, the contribution can be re-expressed as a derivative with respect to the interpo- 
lating parameter An: 



Am(_i,_i) = 2/j2 



dAn 



1 / rn? + An ■ m^ 

1 + An ■ Ao M 1 + An ■ Ao 



(37) 

An=0 



Now, it should become clear how the quenched infrared singularity arises. While the contribu- 
tions from both the pion pole and the shifted pole have the same functional form as the one 
denoted as Xi, in the fully quenched limit, the shifted pole returns back to the unshifted posi- 
tion and gives rise to a derivative contribution, which is not of the form Xi. In fact, this yields 
the precise statement of Bernard-Golterman's third theorem for baryons in the fully quenched 
limit. Expanding the derivative in Eq. (^) explicitly, one finds that 



d 



dAn 



1 jm'^ + An ■ rriQ 

■ Xi 



1 + An-Ao \ 1 + An-Ao J J^^^^ 



1 Q r 

= -(--ml-m Ao-m^)=l2(m'^). (38) 

Note the appearance of the term linear in m. This is the anticipated contribution of non-analytic 
infared singularity in quenched QCD. To summarize, we conclude that the mass correction in 
the fully quenched limit is reproduced correctly and is given by 

Am = / (^AT^ + 1 J2) . X,(m') + 2j2 • Um'). (39) 
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3.5 Large- A^c Suppression of Am(_i _i) 

As mentioned before, the hairpin diagram is OZI suppressed. As a result, the hairpin parame- 
ters, Aq and itIq, are of order l/Nc- In hght of this piece of new information, we will reanalyze 
the 1/iVc and An dependence of A?7i(_i _i), the term which originates from the hairpin propa- 
gators. 

To make the Nc dependences explicit, rescale the hairpin parameters as 



rrin 



mi/N, 



Ao = Ao/N, 



(40) 



where rn^ and Aq have smooth non-trivial large N^. limit. Then the mass correction becomes 



Am,_.-i)=/J=-(^)- 



1 + (An/AfJ^o, 
This can be expand in a Taylor series of An/Nc, 



fm^ + {An/Nc)rhl^ 
'[l + {An/NMo , 



Mm^) . (41) 



2 J^ ^ (An/N,)'' 
Am(_i,_i) = 2/ — 5:-^ / "^ 



An 



2t2 



fc=l 






l2(m' 



k\ d{An/Nc)^ 

An- , n, 
H ^z{rrr) + ■ 



'm^ + {An/Nc)rnl'^ 
1 + {An/Nr)Ao" V 1 + (An/ATjAo , 



^^i 



J (An/Afc)=0 

(42) 



where X2_3 are defined as 



Mm" 



d 



d{An/Nc) 
127r/2 \2 



1 



1 + (An/ATjAo 
5 
2 



'[ 1 + {An/N,)Ao , 



■ TTin ■ m — - ■ Aq ■ m^ 



(An/7Vc)=0 



(43) 



XJm' 



rf2 



1 + {An/N,)Ao 



^■Ii 



diAn/N^y 



1 /3 _4 1 15. 2 

127r/2 V4 ° m 2 



'm^ + (Ara/A^J^o' 
^ 1 + {An/Nr)Ao ^ 



10 ■ /Aig 



771 



An ■ m^ 



J {An/N^)=0 



(44) 



From Eq. (^) one can see clearly that the leading t]' correction is of order g"^ / pNc ~ l/A'c > 
viz. suppressed by one higher order in l/N^ than we have naively expected. Moreover, all 
the higher order corrections are in positive powers of An, which vanish identically in the fully 
quenched limit. The real significance of this result is that, in the leading order of an l/N^ 
expansion which keeps only the first term of Eq. (42), Am(_i^_i) is of the form of X2(m^), a 



form which leads to new chiral divergences such as nn}J'^^ even when the theory is far from 
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the fully quenched limit An ^0! In other word, through taking the large A''^ limit, one can 
bypass the Bernard-Golterman's third theorem, and have quenched chiral singularities even if 
the theory is only partially quenched. Note that our result is not in contradiction with the 
Bernard-Golterman's third theorem: the analysis of Ref. is for A^c = 3. In fact, the third 
theorem holds for any fixed value of Nc but fails only when Nc is taken strictly to infinity. It 
is important to understand that these quenched infared singularities appear when the shifted 
pole coincide with the pion pole. Since the shift of the pole is proportional to An/Nc, these new 
chiral divergences appear both at the fully quenched limit and also at large (but not infinite) 
Nr hmit. 



'c 



We note in passing that the higher order terms in the Taylor expansion Eq. ([42D are more 
singular in the chiral limit. In the chiral limit, while Xi itself goes like m^, X2 and X3 behave 
like m and 1/m respectively. It is easy to verify that each term will be more singular by a 
factor of 1/m?. If we attempt to take the chiral limit in expansion Eq. (|42D before taking the 
large Nc limit, all the terms (except X2) will diverge, and the series will be a sum of infinitely 
many infinities. However, all these divergences are not physical, as the whole series can be 
exactly resummed into the closed form expression Eq. (0). The infinities just assert that 
the true expansion parameter in the Taylor series is {An/2Nc){rn^/m?), which diverges in the 
chiral limit. As a result, the Taylor expansion has behaved badly. In other words, the chiral 
divergences oi Tk{k > 2) are merely an artifact of applying an expansion outside its radius of 
convergence (zero in this case), while all physically relevant quantities can be obtained from 
Eq, (0), which has indeed a well-defined chiral limit. 

This above discussion poses an important implication on whether we can really see the 
m ~ my^ quenched chiral singularities from lattice data. The X2 term dominates the Taylor 
expansion only when rr? > {/\n/2Nc){ffi^) = Anm'^/2, and Am(_i_„i) is linear in m. With 
small pion mass m? < Anmg/2, however, Ar7i(_i „i) will deviate from linearity in m, as the 
higher X are no longer negligible. On the other hand, Ar7i(_i _i) is usually dominated by 
Am(+i_o) except at very low pion masses, i.e., w? < SJ'^ml/N^. So we can observe Am with a 
linear dependence on m only over the window 

Anml/2 <m^ < SJ^m^/A^'f , (45) 

and such a window exists only if 

An < IQJ^Nl (46) 

This condition is independent of the value of the parameter ttiq, and is marginally satisfied for 
An = 1, A^c = 3 and J = |. While this is only an order of magnitude estimate, it suggests that 
it would be difficult to see the linear dependence; we expect that of the Bernard-Golterman's 
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-Nl + -nN^ + (- - -)J^ ■ Xi(m2) + J^ . L — - . J^ ' o 



third theorem is bypassed only when A^^ ^ oo, but in that hmit the right hand side of Eq. 
goes to zero, and the window just disappears. Physically it means that, when A^^ — > cxo, the 
linear term is overwhelmed by the leading order term Am(+i q), which is larger than the linear 
term by a factor of A^^. We conclude by restating that, while the Bernard-Golterman's third 
theorem is bypassed in the large Nc limit, it would probably be difficult to observe in real lattice 
data. 



3.6 Unquenched QCD Limit — /c = 

Let us now consider the conventional QCD with no ghost quarks by taking the limit fc — > 0, 
viz. An = n. We continue working on theory with two flavors only. Thus, Eq. (|30|) becomes 

^ '" " A '2 n' J ' ' 1 + n-Ao \ 1 + n ■ Aq 

(47) 

Note that the term proportional to Xi(r7i^) is exactly the loop correction with one-pion exchange 

in the standard n = 2 xPT. The interpretation of the shifted pole is in the second term is that 

it is the reinstatement of the contribution of rj' meson. The shift of the pole just reflects the fact 

that, in the real world, the t]' mass is shifted from the pion mass due to the nonperturbative 

resummation of necklace diagrams. 

To see this clearer, let us rearrange Eq. (|30D in the following way: 

Am = /-( ^N^ + ^nN, + i^-^)j' y Mm') + g%s ■ Mml>), (48) 

in which 

al'BB = 2 3' g^/Zn, mj, = {m' + n- ml)/Z, where Z = {l + n- Aq). (49) 

These represent coupling, mass and wave function renormalizations of the t]' meson. Note that 
the expression for nir,' and Z agrees with the standard xPT results in Eq. (|]). 



The second term in Eq. ( ^Sf ) represents the contribution of T^'-meson to the chiral one-loop 



correction to the baryon mass spectrum. We thus conclude that the unquenched limit k = 
does not reduce to the xPT in which only octets of Goldstone mesons are retained. Instead 
we have reproduced the xPT without integrating out the rj' meson. This alerts us that we have 
to be careful when applying the Bernard-Golterman's first theorem, which states that in the 
subsector where aU valence quarks are unquenched, the SU(n|fc) theory is completely equivalent 
to a normal, completely unquenched SU(n — k) theory. We found that this is indeed true, but 
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unlike the standard xPT, this "completely unquenched SU (n — k) theory" contains an rj' meson 
which may (and does) contribute to chiral loop corrections. As a result, naive comparisons 
between results of (P)QxPT and their counterparts in standard xPT may be problematic or 
misleading. One should instead compare results of (P)QxPT with their counterparts in xPT 
with rj' loops. 

In passing we observe that it is impossible to compare results of standard xPT with their 
counterparts in (P)QxPT with the rj' integrated out. In (P)QxPT, the rj' propagator always 
has a single pole at m^ except for the unquenched limit k = 0. Consequently, there is no 
?7'-impossible to integrate out the r]' meson. So the scheme suggested above, namely comparing 
results of (P)QxPT with their counterparts in xPT with r]' loops, is the only option which is 
theoretically well-justified. 



4 Discussion 

4.1 1/Nc Corrections 

So far, we have kept only contributions that are leading order in l/N^ expansion. In this 
section we will discuss several sources of higher-order l/N^ corrections we have ignored so far. 
First of all, there are corrections to the axial current couplings to baryons through multi-quark 
operator. It has been shown [^ that the leading order correction to pion coupling is through 
the two-quark operator 



1 n' ^'^ 

- I {d^n'^) Y: {qliS' ® l)q,) {4(1 ® T»)g,) , (50) 

where g' is an undetermined coupling constant of order 0{N^). This will lead to new single- 
pion loop corrections: a gg' / P term at the order of 0{N^^), and a g'"^ / P term at the order 
of 0{N~^). Since these terms are sub-leading with respect to the leading pion contribution 
~ 0{Nc)., they can be safely ignored in the large Nc limit. What cannot be ignored is the 
rj' counterpart. By planar symmetry, rj' can also couple to a baryon through the following 
two-quark operator: 

1 n' ^'^ n' ^^ 

which is of the same form as the leading order rj' coupling. In a consistent 1/Nc expansion, 
one has to keep the effect of g\ which contribute in leading order {0{N~^)) to the quenched 
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Am(_i,_i) = {g + g'f3^ — {l,{M^) - l,{m')) ^ O [N^'/^n-') . (52) 



infared divergences Am(_i _i)Q : 

2 

An' 

This does not change the form of the new chiral singularities, but its coefficient is no longer 
related to the standard xPT correction (Am(i^o))- 

Another possible source of correction is the inclusion of the hairpin coupling h, which is again 
of order l/Nc- We have not calculated Am in the presence of a non-zero hairpin coupling, which 
should be straightforward but adds significantly more complications. However, the functional 
form of Am is not affected by the inclusion of h, as the contributing Feynman diagrams are 
still of the same form as those shown in Fig. 1. Hence, we expect that only the overall factor 
g"^ will be modified: 

g^^(g + xhy. (53) 

Unfortunately, the variable x may be different for different contributions for Am, i.e., x for 
Ar7i(+i 0) may be different from that of Ar7i(_i _i). However, since h is l/N^ suppressed with 
respect to g, the change will only be higher order effects in a 1/Nc expansion. In any case, 
our analysis of the fully quenched and the unquenched limit depends only on the functional 
form taken by Am. The coupling constants are irrelevant for this analysis, and we expect the 
reported results remain valid even with a non-zero hairpin coupling, even though the details 
have to be checked. 

4.2 Conclusion 

In this paper, we have studied baryons in partially quenched large-A^c QCD with particular 
emphasis to the interpolation between fully unquenched and fully quenched limits. In the large 
Nc limit of partially quenched QCD, we have calculated Am, the chiral one-loop correction to 
baryon masses. The main results are: 

• The leading order contribution to Am is of order 0{Nc), in agreement of the large Nc 
counting rules. This contribution does not depend on the value of An, which survives both 
in the fully quenched and the unquenched limits. In other words, the leading term of Am is 
independent of number of ghost quarks introduced. This provides circumstantial evidence that 
quenching correction to baryon masses is small. 

• For any fixed Nc our expression for Am does satisfy the Bernard-Golterman's third 
theorem and has no quenched infared singularities except An -^ limit. In the large Nc limit, 

® The coupling g' appears also in the other contributions to Am, but the effects are subleading. 
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however, the third theorem is bypassed. We indeed find appearance of the quenched chiral 
corrections, which are unfortunately of order C(A^~^) and l/N^ suppressed with respect to the 
leading order contribution to Am and hence making them difficult to be observed. 

• In the unquenched limit, xPT results are reproduced with rj' loop contributions. 

Unfortunately, it is highly non-trivial to relate our result with that of Ref. |^, which has 
studied QxPT for baryons with n = fc = 3 and A^c = 3. The reason is that it turns out that 
the spin-fiavor structure of large N^ baryons with n = 3 are much more complicated than that 
in n = 2. In particular, the large N^ limit of the ratios between the SU(3) coupling constants 



JF, P, C and Ti in Ref. |g] are ambiguous ||22|. The study and clarification of these issues are 
straightforward, but is beyond the scope of the present paper. Being that the motivation and 
the focus of our study in the present paper are of more theoretical issues than phenomenological 
ones, detailed comparison of our results with those of Ref. § will be relegated in a separate 
paper elsewhere. 

We will end the paper by discussing possible relevance of our studies. Our motivations are 
mainly theoretical and we have tried to disentangle the interplays between chiral dynamics, large 
Nc expansion and quenching. On the other hand, while most of the lattice QCD simulation are 
performed with N^ = 3, there are also simulations with arbitrary Nc, most notably by Teper 
[p9| . So far, these studies have focused on pure Yang-Mills theory and have studied string 
tension and glueball masses. However, new data on large Nc QCD should become available in 
the future. More phenomenologically, there has been a lot of investigation of baryon properties 
in the large Nc limit [^, ^, ^ ^ ^ ^ , with many new interesting results. For example, the 
A-N splitting is expected to be of order l/Nc- While these results are successful in organizing 
the baryon properties, one would like to directly verify these predictions with the A-N splitting 
for different values of Nc. This, however, is only possible with lattice simulations. Our studies 
should be relevant if these simulations are really undertaken in the foreseeable future. In 
conclusion, we are using the l/Nc expansion as a guiding principle to help identifying and 
organizing spectrums and chiral corrections when one works on quenched lattice world, hoping 
that l/Nc expansion is as successful in hadrons phenomenology on the lattice as it has been in 
the real world. 
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